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Abstract Induced emission due to the presence of soft CMB photons slightly increases the two-photon decay rate 
of the 2s level of hydrogen defining the rate of cosmological recombination. This correspondingly changes the 
degree of ionization, the visibility function and the resulting primordial temperature anisotropies and polarization 
of the CMB on the percent level. These changes exceed the precision of the widely used Cmbfast and Camb 
codes by more than one order of magnitude and can be easily taken into account. 
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1. Introduction 

One of the key processes for the formation of the pri- 
mordial temperature fluctuations of the cosmic microwave 
background (CMB) is the recombination of hydrogen, 
which at redshifts z ~ 1100 makes the Universe trans- 
parent for CMB photons. The time and duration of 
recombination directly influence the characteristics of 
the CMB anisotropies. Today, line of sight Boltzmann 
codes l ike Cmbfast llSeliak fc Zaldarriaeal HQOfil and 
Camb ijLewis et all l2000|) are routinely used to calcu- 
late the power spectrum of the primordial temperature 
anisotropies with inclusion of many different physical pro- 
cesses in the early Universe. A precision of the solution 
to ^ 0.1% within the assumed models up to multipoles 
I ^ 3000 is reached. This level of precision is now becoming 
necessary with the advent of space missions like Wmap, 
Planck and ground-based experiments such as Act and 
Spt, which will allow us to measure the CMB tempera- 
ture and polarization anisotropies with unprecedented ac- 
curacy and thereby open a possibility to determine the 
key paramet ers of the Universe with high precision. 

Recently iDubrovich fc GrachevI l)2005l) have included 
the two-photon decays of high levels of neutral hydrogen 
and helium in their calculation of the recombination rates, 
in the case of hydrogen yielding corrections on the level 
of a few perce nt and a sig niflcant acceleration of helium 
recombination. iLeung et al. (2004) have included the soft- 
ening of the matter equation of state due to the transition 
from completely ionized to neutral matter and found that 



the CMB temperature and polarization power spectra are 
affected on the level of some percent at large multipoles. 
As they pointed out, these corrections exceed the level 
of cosmic variance at multipoles / ^ 1000 and therefore 
should be taken into account in high accuracy analysis 
of the CMB data. In this paper we discuss an additional 
physical process that changes the ionization degree of hy- 
drogen in the Universe at any given moment of recombi- 
nation on the level of a few percent. 

It is generally accepted that the rate of recombina- 
tion is mainly controlled by the t wo-photon decay of the 
metastable 2 s level of hydrogen llZeldovich et all Il968l: 
IPeebleslllQeS ISeaeer et al.1ll999l I^OOO"). Here we discuss 
the influence of the simula ted two-photon emission (e.g. 
see lBerestetskii et al.ll97ll p. 229) due to the presence of 
the low frequency photons of the CMB blackbody radia- 
tion. Below we present a simple calculation for the change 
of the two-photon decay rate of hydrogen, A2sis, with red- 
shift and shortly discuss the consequences for the visibility 
function and CMB temperature and -E-mode polarization 
power spectrum. 

2. Induced 2s Is two-photon transition of 
hydrogen during recombination 

Based on the pioneering work of ICoppert-Mave 3 lll981^ 
the rate for the 2s — > Is two-photon transition of hydro- 
gen, assuming no ambient photon fleld, has been calcu- 
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for the total transition rate, which to 0.1% agrees with 
a more simple calcula tion based on the method used by 

zer fc GreensteinI l|l951 ). In this paper we use this sim- 
plified approach and include induced effects in the calcu- 
lation of the hydrogen two-photon decay rate within the 
context of cosmological hydrogen recombination. A simi- 
lar calculation can be used to include the effects of induced 
two-photon decay for helium recombination. 

The total transition rate with no ambient radiation 
field can be given by 



^2sls — 



(1) 



with Aq = 9a^cR/2^^ « 4.3663 s~^, where a is the fine 
structure constant, c is the speed of light and R is the 
Rydberg constant for hydrogen. In equation 0(y) d?/ 
is proportional to the probability of emitting one photon 
at frequency y = v/vq in the range Ay — dv/u^, where 
vq is the frequency of a Lyman-a photon, with energy 
^ 10.2 eV, while the second photon is emitted at y' = 
l — y= [uQ — i']/vf). The factor of 1/2 is required since there 
are two photons and each pair is co unted twice. The func- 
tion M il') is defined in the paper of lSpitzer fc GreensteinI 
lll95lh (cf. Eg. 3), which nowadays can be easily calculated 
numerically. 

In the presence of an ambient radiation field with oc- 
cupation number n{v) the total transition rate is given by 
the expression 



A 



ind 
2sls 



n{u)][l + n{va~v)]Ay, (2) 



where the factors of 1 -I- n account for the effect of Bose- 
bunching. In the context of recombination the radiation 
field is given by a blackbody spectrum for which the occu- 
pation number is n(j/) = i/Je^'^/'^^^ — 1], with the photon 
temperature T = Jn(l + z), where Tq = 2.725 ± 0.001 K 
i Fixsen fc Matheill2002l) . The relation between hv/kjiT = 
X and y is given by 



^ = 1100 

43455 i 
X = — — y « 40 y . 

1 + z 



(3) 



Figure ^ shows the comparison of the integrands of 
and |5J|. Due to induced effect the probability of emitting 
one soft photon {y ^ 0) and at the same time the other 
close to the Lyman-a frequency (y ~ 1) is enhanced. This 
enhancement due to the change of x with time depends on 
the redshift. At sufficiently small redshifts, induced effects 
become negligible. 

One can examine the behavior of 0(y) for y 
1 more closely by using the analytic fit as given by 
iNussbaumer fc Schmut3 l)l984 '). With this, one obtains 



(j){y) = C [w{l -4r'w^) + a w'^'^'' 4r<] 



(4) 



where w = y\l ~ y], C = 46.26, a = 0.88, = 1.53 and 
7 = 0.8 (note that in INussbaumer fc Schmutzl (jl^SJ) the 
normalization constant C was defined differently). With 




Figure 1. Two-photon 2s decay of hydrogen: the solid line 
shows t he two-photon probability distribution, 4>{y), as 
given in lSpitzer fc GreensteinI l)l95l[l assuming no ambient 
radiation field. In contrast to this, the dashed line includes 
the effects of induced emission due to the presence of the 
CMB at a redshift of z = 1500. 

this the integrands of equations and ^ around y ^ 
can be expressed as 



/ind 



0(y)«Cy[l-3.O3yO-8] 
C 



(y) 
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respectively, where here we introduced 
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(5a) 
(5b) 

(6) 



These approximations are accurate within a few percent 
for y < 0.05 in the redshift range z = 1000 - 1500. 

Here it is important that (t>{y) vanishes as (f>{y) ^ y 
for y 0. Since in this limit n(iy) l/y, the product 
(j){y)n{v) is finite. Therefore it is not necessary to intro- 
duce a low frequency cut-off in the integral . 

From ^ one can also deduce the contribution of the 
low frequency part to the total two-photon decay rate. 
Integrating from zero up to ?/ni one finds 



A^2sis « 50.5 y^Ji- 2.16 s'^ 



ind 
2sls 
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(7a) 
(7b) 



Here A indicates that only part of the integrals ^ and Q , 
i.e. for y G [0, ?/,„] have been taken. Due to the symmetry 
of My) around y = 1/2 one can apply the approximations 
(ISJ and Q also to the case y ^ 1 by simply replacing 

If one assumes that photons very close to the center 
of the Lyman-a line cannot escape then this reduces the 
contribution of the two-photon decay to the effective rate 
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Figure 2. Rcdshift dependence of the relative change of 
the total two-photon transition rate of hydrogen, [A2s'j's — 
^2sis]/^2sis with A2sis — 8.2206 s~^, in the context of 
recombination using Eq. ||2Jl. In addition the results ob- 
tained with the simple analytic approximation for (f>{y) 
are shown. 



of recombination, since the trapped photons prevent the 
corresponding recombination. As an example, if at rcd- 
shift z — 1100 photons within 1% (5%) of the central 
frequency of the Lyman-a line cannot escape then the ef- 
fective two-photon decay rate (i.e. the rate used in the 
calculation of the recombination history) is smaller by 
~ 0.06% (~ 1.23%) for the standard calculation of A2sis 
and by ^ 0.33% 2.01%) for the calculation including 
both spontaneous and induced effects. 

As will be shown below, in general the contributions 
due to induced effects are at the level of some percent 
themselves. Hence, this estimate shows that the correc- 
tions to the two-photon decay rate we are discussing here 
can only be considered accurate if the more energetic pho- 
tons from the 2s two-photon decay lying within less than 
~ 0.1 — 1% of the Lyman-a line center are trapped. At 
redshift z ~ 1100 the Doppler width of the Lyman-a 
line is Ai^D ^ 2.3 x lO^^t'o. Our computations of the 
Lyman-a photon escape in the distant wings show that 
there is no significant diffusion back to the line center be- 
yond a few ten to hundred Doppler widths. Therefore the 
aforementioned condition should be easily fulfilled, since 
at z = 1100 a 1% distance from the line center corresponds 
to ~ 435 Doppler widths. However, even if every photon 
within 1% of the line center is unable to escape, this would 
only affect the results obtained below by about ten per- 
cent. One should note that the expected corrections due 
to the Lamb-shift are much smaller. 

In Figure H we present the rcdshift dependence of the 
relative change of the total two-photon transition rate of 
hydrogen in the context of recombination using equation 
A simple fit to this function, which within ^1% ac- 



curacy is applicable in the redshift range 500 < z < 2500, 
can be given by the expression 

^-^-^ = 1.181 X 10-3 ^ + 2.177 X 10-2 

^2sls 

- 1.958 X 10-3 X^ (8) 

with X = (l + ^:)/1100 and where A2sis and Al^^^^ are given 
by equations 1^ and respectively. The main correc- 
tion scales as oc (1 + z)^. The total correction exceeds the 
percent level for z > 700. In addition the results obtained 
with the simple analytic approximation Q for 4i{y) are 
shown. They agree very well with the full numerical result 
but were obtained with much less numerical effort. 

The rate for the inverse process can be found from 
equation Q using the principle of detailed balance. 
In thermodynamic equilibrium the ratio of the popu- 
lation of the Is and 2s levels is given by A^2s/-^is — 
exp{-E2i/kBT), where £'21 = 10.2 eV. Therefore the 
rate for the inverse process is given by Ais2s — 
Ajsis exp(— £'21/^3?"). However, during the period of 
recombination, which is most relevant for the CMB 
anisotropics, the inverse process is negligibly small. 

3. Changes of the ionization fraction, visibility 
function and power spectra 

It is straightforward to include the redshift dependence 
of the 2s two-photon transition rate into the Cmbfast 
code using the approximation (jS]) . One only has to replace 
the standard two-photon decay rate A2sis by A^2sisi^)- 
With this one can calculate the changes in the ionization 
fraction, the visibility function, V{z) = exp(— r) dr/ drj 
l|Sunvaev fc ZeldovichI Il970(l and the temperature and 
E-mode polarizat ion power spectra fo r the Wmap con- 
cordance model (jBennett et al.ll2003|) . Above, r is the 
Thomson optical depth and rj is the conformal time. 

In Figure|31we have presented the redshift dependence 
of the ionization fraction and the corresponding change 
due to the inclusion of induced two-photon emission for 
the Wmap concordance model using a modified version 
of Cmbfast. The ionization fraction is affected by a few 
percent with a maximal relative difference of ~ —1.3% 
a.t z ^ 1050. Because the total two-photon decay rate is 
slightly higher than in the standard calculation, recombi- 
nation occurs a bit faster. 

In Figure 01 we show the results obtained for the visi- 
bility function. At redsh ifts much below the maximum at 
Zdcc = 1089±1 (jBennett et al. 2003.1 the visibility function 
is affected by less than one percent, whereas for z ^ z^cc 
the change reaches 8%. The position of the maximum and 
the width of the visibility function are both affected on a 
level below one percent. 

In Figure|21we give the results obtained for the temper- 
ature and _E-mode polarization power spectra. Again the 
changes due to induced two-photon emission are on the 
level of a few percent, with an increase towards smaller 
scales. The amplitude of the change of the polarization 
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Figures. Ionization fraction, = Nc/Nn, and the rel- 
ative change due to the inclusion of induced two-photon 
emission for the Wmap concordance model. Here Nc and 
-/Vh are the free electron and hydrogen number densities, 
respectively. 
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Figure 4. Visibility function, V{z) — exp(— r) dr/ drj, and 
the relative change due to the inclusion of induced two- 
photon emission for the Wmap concordance model. The 
amplitude of V(z) is normalized, such that the maximum 
is 1. 
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Figure 5. Temperature and i?-mode polarization power 
spectra and their relative change due to the inclusion of 
induced two-photon emission for the Wmap concordance 
model. 

power spectrum is roughly twice that of the temperature 
power spectrum. 

4. Conclusion 

Due to induced two-photon decay of the hydrogen 2s level, 
the rate of recombination is increased on the level of a few 
percent around the maximum of the visibility function. 
This increase results in changes of the ionization fraction, 
the visibility function and the temperature and polariza- 
tion power spectra by a few percent. These changes can be 
easily taken into account for future high accuracy analysis 
of CMB data using the approximation 0. Induced two- 
photon decay would similarly influence the recombination 
of HcII and HcIII, but in that case the effects on the CMB 
power spectra are expected to be extremely small. 
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